Microtubule Length-Regulation by Molecular Motors 
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Length-regulation of microtubules (MTs) is essential for many cellular processes. Molecular mo- 
tors like kinesin 8, which move along MTs and also act as depolymerases, are known as key players 
in MT dynamics. However, the regulatory mechanisms of length control remain elusive. Here, we 
investigate a stochastic model accounting for the interplay between polymerization kinetics and 
motor-induced depolymerization. We determine the dependence of MT length and variance on rate 
constants and motor concentration. Moreover, our analyses reveal how collective phenomena lead 
to a well-defined MT length. 
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During the lifespan of an eukaryotic cell microtubules 
(MTs) perform highly dynamic tasks. For instance, 
during mitosis, they form the mitotic spindle, which 
searches, captures, and separates the double set of chro- 
mosomes [1]. To achieve such complex dynamic behav- 
ior there need to be molecular mechanisms which allow 
a dynamic control of MT length. There is much evi- 
dence that these mechanisms rely on an intricate inter- 
play of GTP hydrolysis [2], mechanical forces [3], and 
regulatory proteins [4]. In particular, the role of the 
molecular motor families kincsin-5 and kinesin-8 has been 
investigated: Several in vivo experiments showed that 
both, the presence and the concentration of such pro- 
teins, strongly affect the functionality of the mitotic spin- 
dle [5, 6]. This is supported by in vitro experiments which 
specifically studied the molecular mechanisms of interac- 
tions between motor proteins and microtubules [7-13]. 
In general, it is accepted that kinesin-8 hampers MT 
growth. In particular, it was found that the plus-end 
directed motor kinesin-8 of budding yeast, Kip3p, de- 
polymerizes MTs at the tip. To gain a deeper under- 
standing for the molecular mechanisms underlying these 
depolymerization dynamics Varga et al. [7, 8] studied the 
interaction of Kip3p with stabilized MTs not exhibiting 
dynamic instability [2, 14]. The key result of these exper- 
iments is that depolymerization is length-dependent, i.e., 
longer MTs depolymerize faster than shorter ones. One 
main determinant of the observed length-dependence are 
molecular traffic jams which can successfully be described 
by driven diffusive processes [15]. These findings suggest, 
that length-dependent depolymerization in combination 
with polymerization allows a cell to regulate the length of 
MTs [7, 8] . There are by now several theoretical studies 
addressing length-regulation ranging from MTs [16, 17], 
over actin filaments [18] to fungi [19] and flagcllae [20]. 

In this Letter, we study the combined influence of 
spontaneous MT polymerization and motor induced de- 
polymerization. In our model we neglect MT dynamics 
at the minus end as there the dynamics rates are much 
smaller than the ones at the plus end [1]. Furthermore, 
under physiological conditions often the minus end dy- 
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FIG. 1. Illustration of the model. Motors attach to and de- 
tach from the MT lattice at rates cu a = cw a and cud, respec- 
tively. On the lattice particles hop to the right at rate v 
provided that the next site is empty. At the right bound- 
ary, the MT plus-end, particles remove the last lattice site 
at rate 8 and the MT lattice grows at rate 7. The resulting 
antenna-like density profile p-(x) is sketched in light gray. 



namics are completely suppressed due to capping pro- 
teins [21]. We build on a recently validated quantita- 
tive model for MT depolymerization [8, 15], and extend 
it by introducing polymerization dynamics at the fast- 
growing plus-end [1] . This accounts for MT growth me- 
diated by spontaneous [2] or enzymatically catalyzed [22] 
attachment of tubulin hctcrodimers to the tip. This ap- 
proach enables us to study the basic principles under- 
lying length-regulation which is achieved by the antag- 
onism between length-dependent depolymerization and 
spontaneous polymerization dynamics. We predict quan- 
titative criteria for the parameter regime where regula- 
tion is feasible. In addition, we calculate both the mean 
length and the corresponding standard deviation, and 
thereby determine the accuracy at which regulation is 
achieved. To describe the MT dynamics we employ a 
driven diffusive lattice gas model [23, 24] as illustrated 
in Fig. 1. Since MT protofilaments serve as independent 
tracks for the motors [25, 26], a MT can effectively be 
described by a one-dimensional lattice of dynamic size 
L(t). The size of a tubulin heterodimer sets the basic 
length scale of the lattice. The state of each site, i, is 
described by its occupation number, n, € {0, 1}, where 
n, = and n, = 1 signify an empty and occupied site, 
respectively. On the MT lattice the dynamics follow the 
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totally asymmetric simple exclusion process with Lang- 
muir kinetics [24] : Motors can attach to and detach from 
the MT at rates uj a = cui a and w^, respectively, where c 
is the motor concentration in the surrounding fluid; the 
binding constant is defined as K := ujjujd- On the lat- 
tice, particles move to right at rate v provided that the 
next site is empty; v = 1 sets the basic time scale. The 
combined effect of motor attachment in proximity of the 
minus-end and subsequent movement towards the plus- 
end leads to an accumulation of motors, which finally 
results in an antenna-like steady state profile [8, 15, 27] 
as illustrated in Fig. 1. At a certain distance from the 
minus-end the density profiles saturate to the equilibrium 
Langmuir density p^n = K/{K + 1) [28]. The resulting 
accumulated density profiles in vicinity of the minus-end, 
p~(x), can be described by Lambert- functions [15, 24]. 
Moving further towards the right boundary (MT minus- 
end) , the density profile is determined by the interplay of 
motor current and the boundary conditions at the plus- 
end. This entails a rich variety of collective phenomena 
and leads to nontrivial density profiles [29, 30]. In the 
present study, the right boundary is dynamic. Motivated 
by the recent studies on kinesin-8 [6-12], we consider the 
following scenario: When a motor arrives at the MT tip, 
it detaches by removing the last MT-site at rate S [15]. In 
addition, subsuming the effects of spontaneous and enzy- 
matic polymerization, the MT is assumed to polymerize 
through the attachment of single tubulin hcterodimers at 
an effective rate 7. These boundary conditions lead to a 
dynamic MT length which is determined by the combined 
effect of the particle current onto the last site, polymer- 
ization, and depolymerization rates. 

The dynamic length of the MT, L(t), is determined by 
the particle density at the MT plus-end p+(L), 

d t L(t) = -Sp+{L)+r (1) 

This equation defines a critical density p c + ="//$, at which 
the MT length is in a steady state, dtL = 0. For tip den- 
sities smaller or larger than p c + the MT grows or shrinks, 
respectively. As the tip density is fed by the motor cur- 
rent towards the tip, it depends on the accumulated mo- 
tor density in bulk p~{x). This suggests the following 
mechanism for MT length-regulation: On short MTs, 
the accumulated motor density is low, and therefore also 
the tip-density p + (L). As long as p+{L) < the MT 
grows. In contrast, for longer MTs higher accumulated 
motor and tip densities are reached which eventually re- 
sult in MT depolymerization once p+{L) > p c + . However, 
this mechanism is only expected to work if the tip is not 
growing too fast: Above a critical polymerization rate 
the particle current feeding the tip density can no longer 
follow the advancing tip. 

To quantify these heuristic arguments and determine 
the precise conditions under which length-regulation is 
feasible and which length is adjusted, the tip density has 



to be determined. This requires to analyze the intri- 
cate interplay between molecular crowding due to high 
motor density [23, 24] and transport bottlenecks at the 
plus-end [15, 31]. In addition, this boundary is highly 
dynamic, and calculations of the tip density are more 
intricate than for standard driven diffusive models for 
which the size of the lattice is constant [29, 30, 32]. 

To make further progress, we first consider a simplified 
model (SM) where we disregard spatial variations of the 
density profile. In detail, we assume a constant density 
p- that serves as a particle reservoir at the left boundary, 
neglect attachment and detachment kinetics, but leave 
the dynamics at the plus-end unchanged, see Fig. 2(a). 
This allows us to focus on the dynamics at the plus-end 
and to unravel how they depend on the reservoir density 
p_. Since, we find that the density profiles adapt adia- 
batically to a dynamic lattice size [33] , the results for the 
full model can be inferred upon replacing p_ by p~{x). 
As the length of the lattice is dynamic, we perform our 
calculations in a comoving frame fixed to the right bound- 
ary. In this frame, a polymerization event corresponds to 
the simultaneous movement of all particles on the lattice 
to the minus-end by one unit, while depolymerization re- 
sults in an instantaneous shift to the right. Thus, in a 
mean-field approximation [(niUj) = (rii)(nj) = piPj] the 
particle current in bulk is given by, 

J(Pb,P+) = Pb(l - Pb) - jPb + Sp + p h , (2) 

where pb is the motor density in bulk. The first term de- 
scribes the hopping processes, the second and third term 
account for simultaneous movement of all particles due to 
polymerization and depolymerization, respectively. Im- 
portantly, the bulk current explicitly depends on the tip 
density and thereby on the right boundary. 

To determine the phase behavior we employ the Ex- 
tremal Current Principle (ECP) [34-36] relying on two 
velocities: The collective velocity v co n(p) = d p J deter- 
mines the direction in which a local density perturba- 
tion spreads. Thereby, one is able to determine whether 
a certain bulk density is stable against perturbations, 
i.e. for a density p stable at the left (right) boundary 
v co ii(p) > (fcoii(p) < 0) holds. The boundary con- 
ditions result in densities at the plus and the minus- 
end, respectively, whose stabilities can now be tested em- 
ploying w co ii. If these densities are stable against small 
perturbations, we call them p lclt and ( o ri s ht as they are 
given by the system's left and right boundary, respec- 
tively. If cither one or both of these boundary densi- 
ties are not stable, perturbations change these densities 
and p leit and p rlght are given by the first stable density 
which is determined by v co n(p) = 0. The shock veloc- 
ity v shock (p lctt , P ^ ht ) = (J(p lctt )-J( P ^ ht ))/(p lctt -p^ ht ) 
determines the direction in which a virtual domain wall 
between the densities at the left and the right, p left and 
p nght , moves and thereby which of these densities is real- 
ized in bulk. In more detail, for u s hock > the left den- 
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sity, p lcft , dictates the bulk density, while for -u s h oc k < 
the right density, p rlght , is realized. In our model particles 
are transported to the right and therefore jams spread 
from right to left. Hence, the virtual domain wall arises 
at the right boundary and the tip densities p!f; ft and p" ght 
determine v s hock; see Supporting Material. 

Due to particle conservation at the plus-end of the MT, 
d t p+ = J(pb, P+) ~ the stationary value of the bulk 
and tip density are related through 



p+5 = J(p b ,p+) 



(3) 



This implies that the values of these densities as well as 
the nature of the ensuing non-equilibrium steady state 
are strongly interlinked. In general, one expects three 
phases [34, 37]: the steady state may either be domi- 
nated by the motor densities at the plus-end (EX-phase) 
and the minus-end (IN-phase), respectively, or by the 
transport capacity (maximal current) of the lattice it- 
self (MC-phase). We first consider the IN-phase where 
Pk N = p_ holds, and Eq. (3) leads to the tip density 



pf(p-) = P-(l-p-- 1 )/{S(l-p_)} 



(4) 



This solution is stable against perturbations only if the 
collective velocity v co i\(p-) = d p J(p, p™ (p))\ p=p _ is posi- 
tive, which holds for reservoir densities smaller than the 
bulk density in the MC-phase p^ c = l—^Fj. If the reser- 
voir density exceeds this value, the ECP implies that 
the tip density becomes constant and independent of the 
reservoir density p+ c = (l — ^/) 2 /5. For the EX-phase, 
the right boundary determines the bulk density p^ x = 
pf x , and Eq. (3) leads to p+ x = 1-7/(1 - S). According 
to the ECP, this solution is stable if the corresponding 
collective velocity w C oii(/0+ x ) = d p J(p, /0+ x )| p=p ex is neg- 
ative. Since in the relevant parameter regime u co ii < 
is always fulfilled, the density p£ x is always stable and 



right 
P+ : 



= pi? x holds. 



In summary, we have found the following results for 
the densities at the left and right boundary of the MT: 



p!f = Min[p I + N ,P+ C ] 



4 sht = Pf ■ 



(5) 



With these expressions at hand, we can now map out 
the phase diagram upon evaluating the shock velocity 
v s hock{p+ ft ,P+ sht ), cf. Fig. 2. The IN-phase is deter- 
mined by 7 < (1 — p-) 2 and S > p_. Importantly, it is 
the only phase in which the tip density is a function of p- ; 
see Eq. (4). As p_ corresponds to the spatially varying 
density profile p~(x) in the full model, length-regulation 
is feasible in this range of parameters. In contrast, in 
the EX-phase [7 < (1 - 8) 2 and 5 < p_] and the MC- 
phase [7 > (1 — S) 2 and 7 > (1 — p-) 2 ] neither the tip 
nor the bulk densities depend on p_. To confirm these 
and the following analytic results, we performed exten- 
sive stochastic simulations employing the Gillespie algo- 
rithm [38]. For both, the simplified and the full model 
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FIG. 2. (a) Illustration of the SM. (b) Phase diagram as a 
function of the depolymerization and polymerization rates, S 
and 7, respectively. The gray shaded area indicates regions in 
phase space in which regulation is possible in the full model. 
The gray area indicates regions in phase space where the MT 
shrinks in the SM. In the full model, for p_ = pLa length 
regulation is only possible in the gray area as detailedly ex- 
plained in the main text, (c) Comparison of simulation data 
with analytical results for the tip density p+. 



discussed in the following, calculations are in excellent 
agreement with simulations, cf. Fig. 2(c) and Fig. 3. 

Moreover, upon combining the results for the tip den- 
sities in the various phases with the critical density 
p c + = 7/(5, we are able to calculate the critical growth 
rate j c , at which the MT length becomes stationary: 



lc(P-) 




EX-phase , 
IN-phase , 
MC-phase . 



(0) 



For 7 > 7 C the lattice grows to infinity, while it shrinks 
indefinitely for 7 < j c . 

Up to now, the discussion was restricted to a simpli- 
fied system, and we have learned how a constant reser- 
voir density translates into the tip density and in which 
parameter regimes the MT grows and shrinks, respec- 
tively. In the following we transfer the so far obtained 
results to the full spatial model, in which the reservoir 
density is replaced by the density profile: p_ — s- p-(x). 
This implies that also the tip density becomes length- 
dependent, p + —> p+{L) in the IN-phase, see Eq. (4). Let 
us first consider how these spatial density profiles affect 
the critical growth rate, y c , in the full model, and thereby 
derive a condition for the parameter regime where length- 
regulation is feasible: Growth is unbounded only if the 
highest accumulated density p—(L) does not result in 
strong enough depolymerization dynamics to overcome 
MT growth due to polymerization. Recall that the accu- 
mulated density profile increases from left to right until 
it saturates to the Langmuir density, pL a - Thus, growth 
is unbounded if 7 > 7 c (pLa)- In contrast, in the de- 
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FIG. 3. (a) Kymograph (upper left): Molecular motors (green shaded traces) accumulate along the lattice (gray) resulting in 
a steady MT length (dashed). Corresponding potential U(L) (bottom left) and length distribution V(L) (right) for K = 1.5. 
(b) Analytical (lines, Eq. (8)) and numerical results (symbols) for the typical MT length L* are compared, (c) Standard 
deviation of the MT length a in units of its typical length L* for the same values for 7 as in (b). Inset: V(L) for 7=1/4 shows 
an exponential tail. [Parameters: 5 = 0.5 7 C = 0.25 (a-c), o; £ j = 2xl0 -4 (a), and lj^^IxIO -3 (b), (c)]. 



polymerizing regime, 7 <7c(/?La); the MT shortens until 
the tip enters the antenna profile. Within this regime, 
the accumulated density and thereby the tip density de- 
crease with every depolymerization event until the MT 
length reaches a stable fixed point L* , at which growth 
and shrinkage balance each other. As the correspond- 
ing restoring force is conservative, the length regulation 
dynamics can be described by a potential U. It follows 
from —8lU = —8p+(L)+j and leads to an adjusted 
length fluctuating around the mean, as observed in the 
MT dynamics; see Fig. 3(a). 

To calculate the adjusted MT length L*, the full spa- 
tial density profile p- (x) as obtained from mean-field the- 
ory [24], and the stochastic growth and shrinkage have 
to be considered. They can be combined in an effective 
master equation, where the degrees of freedom from the 
occupation numbers, rij, are adiabatically eliminated: 

d t V(L)= [(E+-l)*p + tL) + (E--l) 7 ] V(L). (7) 

Here, are step operators which increase or decrease 
the lattice length; p+{L) is the density at the tip de- 
pending on L. In the IN-phase, in which regulation is 
feasible, P+ (L) = p_(L)(l - P -{L) - j)/[8(l - p_(L))\ 
holds, where p- (x) is the spatial density profile given 
by Lambert- W functions [15, 24]. We solve the Master- 
equation approximately using the van Kampen system 
size expansion [39]: The deterministic dynamics £(t) is 
separated from the fluctuations £ employing the ansatz 
L = fl£(t) + \/f2£. As expansion parameter we consider 
fi = 1 /uj a because the typical length scale of the accumu- 
lated density profile which triggers length-regulation is 
given by l/w a . Additionally, time has to be rescaled ac- 
cording to r = uj a t since the equilibration time also scales 
with this length scale. An expansion of Eq. (7) in terms 
of 1/y/n yields the mean MT length 

L* = (l- U In I ^> v™~i |) . 

(8) 

As can be inferred from Fig. 3(b), this result is in excel- 
lent agreement with numerical data. We observe that the 



stationary MT length is independent of 8, and a mono- 
tonically decreasing function of the binding constant K. 
The latter behavior reflects the increase of the slope of 
the antenna profile with larger K implying that the den- 
sity at which regulation arises is reached for shorter MTs. 
The van Kampen approximation also gives the variance, 



uj a -14- ^/T - 47" +2 7 (1 + if)' w 

For small values of 7, the standard deviation a is be- 
low 10% of the filament length in a range of approxi- 
mately 1 pm— 20 pm. The variance actually diverges with 
7— > PLa(l — PLa) for K < 1, while for K > 1 regulation 
remains possible for 7 = pLa(l — Pha)- In this regime, 
the MT length distribution V(L) develops an exponential 
tail. This tail cannot be described by the van Kampen 
expansion, which explains the deviations between the nu- 
merical and the analytical results in Fig 3(c). 

In this Letter, we investigated how motor-induced de- 
polymerization in combination with spontaneous poly- 
merization can result in length-regulation of biological 
filaments. We found a broad parameter regime in which 
length-regulation is feasible, due to collective phenom- 
ena of molecular motors which also act as depolymerases. 
Even though the regime where length-regulation is pos- 
sible depends on the depolymerization rate, the adjusted 
filament length is independent of the depolymerization 
rate 8, because of microscopic traffic jams forming at 
the tip. Our model provides a proof of principle that 
spatial dependences in the growth and shrinkage rates 
of filaments, which arise from motor transport in this 
case, can result in a well-defined filament length. It may 
serve as a basis for mechanistically more detailed analyses 
which account for multiple protein species [40] , dynamic 
instability [41, 42], internal states of MTs or motors [43], 
assemblies of MTs [44], or the abundance of molecules 
in the cell [45]. We expect, however, that the main idea 
- feedback between polymerization dynamics and collec- 
tive motor dynamics - remains the core mechanism. 
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In this supporting information, wc give an intuitive argument for the extremal current principle to provide some 
background which facilitates understanding. Furthermore, we explicitely show, why the density at the right boundary 
is always stable in the relevant parameter regime. 



PHASE BEHAVIOR 



An Intuitive Argument for the Extremal Current Principle 



To decide which phase is realized in bulk the extremal current principle (ECP) can be used [1-4]. Employing the 
shock velocity, w s hock one can derive whether the right or the left boundary determines the phase behavior. This 
means that the phase transition lines between IN- and EX-phasc or MC- and EX phase can be calculated. Note, that 
the second order transition between IN- and MC-phase is due to the collective velocity (as discussed in the main text) 
and therefore the description employed here does not apply. In the following, wc provide an intuitive argument based 
on traffic jams on the lattice, to decide which phase is realized. As it will turn out, this leads to the same result as 
obtained from the shock velocity. An additional merit of the heuristic arguments is that they provide insights into 
why it is necessary to consider the tip densities instead of the bulk densities in the shock velocity, i>shock(p+ ft > P+ sht )- 
For simplicity, let us assume p_ < p ma x, i-e. the maximal current never determines the transport in bulk. Thus, 
we focus on the IN- and EX-phase and the corresponding phase transition. The regime p_ > p max and the phase 
transition from MC- to EX-phase can be analyzed analogously simply by replacing p_ by /3 max - For p_ < /5 max , there 
are two possible scenarios for the bulk densities as sketched in Fig. 4: In the IN-phase (blue line), the density at 
the left boundary is given by p_ and the tip density has a distinct value, p± , which lies either above or below p_. 
In contrast, in the EX-phase (red line), both the bulk and the tip density are given by the same value p+ x . Before 
turning to the question which phase is realized depending on the parameters, we first show that either both possible 
tip densities lie above p_ (indicated by a '>' subscript), /0 I ) N > and p^y, or below /?_ (indicated by a '<' subscript), 

pf, 

< and p^^. Upon employing Eq. (4) from the main text for the tip density, p- ^ can be expressed as, 

6 1 -p_) 
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FIG. 4. Illustration of the right boundary determining the phase behavior for bulk densities, p_. Density profiles corresponding 
to the IN-phase are shown in blue while the ones corresponding to the EX-phase are red. The tip densities can either both 
be smaller or larger than the bulk density, p™<, P+ X < and p+^>, p+^>, respectively. Both scenarios are depicted here. If the 
dynamics are determined by the entering current as in (A), the tip density, which would arise in the EX-phase, p+ x , serves as 
a virtual bottleneck. Therefore, the tip density corresponding to the maximal exiting current is realized. In other words, the 
IN-phase arises for p™< > p+ x < or p™> < p+ x > (A), while the EX-phase emerges for p™< < p+ x < or p™> > p+ x > (B). 
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Rearraging yields, 

9- S l - jh s = pI x - (ID 

With these results at hand we can now decide which phase is realized in the system. To this end we need to compare 
the current from the left [main text Eq. (1), J loft = J(p~, and the current from the right [J ri s ht — J(p^ x , p+ x )]. 
If p 1 + < > or p I | N > < p+^> holds, the particle current from the left, which depends on p™ [main text Eq. (1), 

J{p~i p l +)\, is always smaller than the one from the right, which depends on p+ x [J(p5 , p5 )]• Therefore, no traffic 
jams arise at the right boundary. Therefore the density profile resulting from the left boundary is not disturbed and 
the IN phase is present, see Fig. 4A. In contrast, for p I + < < /?+^; or p 1 + > > p+ >, the tip density determined by the 
exiting current, p+ x , acts as a bottleneck and a traffic jam results at the tip. As this traffic jam becomes macroscopic, 
i.e. it spreads back into the system it results in a bulk density given by the tip density p5 x ; the system is in the 
EX phase (Fig 4B). Analogously, one can analyze the phase transition between the MC- and EX- phase for Langmuir 
densities larger than p max . 



STABILITY OF THE RIGHT BOUNDARY DENSITY p+ ight 

In analogy to the treatment of the left boundary in the main text, we here analyze the collective velocity at the 
right boundary to decide whether p^ x is stable or unstable. The collective velocity is given by 

v c (pf) = d p J{ PlP ^)\ p=pT = 1 - 2p* x - 7 + 8p* x . (12) 

Thus the p+ x is stable (unstable) if pf x > (p+ x < ^) holds. Equating this condition with p^ x = 1 - 7/(1 - 8), 
one calculates in which parameter regime the density at the right boundary is stable. Thereby one arrives at the 
condition, 

7<(1-<S) 2 . (13) 

This is identical to the boundary between the EX and the MC phase. Therefore, the parameter regime where the 
p^ x is unstable is not important as the p+ then determines the system, 
merlin.mbs apsrcv4-l.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked 
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